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It is found that when gravitational effects are taken into account, the two historical independent tradi- 
tional hallmarks of ideal superconductors—perfect conductivity (a) and perfect diamagnetism charac- 
terized by zero permeability (u)—are actually related to each other in terms of a theoretical o-u model 
developed in this paper. Our result suggests a possible importance of gravitational effects in understand- 


ing superconductivity. 


I. INTRODUCTION 


Important space and other possible gravitational ex- 
periments are crucially dependent on the use of a super- 
conductor shielding system.'~? Experimental efforts to 
realize idealized near-zero magnetic fields inside super- 
conducting shields, such as by using nested and mechani- 
cal expandable superconducting shields, have yielded 
remnant internal magnetic fields of the order of 107 '? t.4 
The attenuation of magnetic fields by superconductors, 
which is a subject of fundamental importance that as yet 
has not been satisfactorily explained, has recently become 
critically important to understanding the interaction be- 
tween gravitation and superconductivity. We have re- 
cently demonstrated that a residual magnetic field, given 
by Eq. (31) in Ref. 5 as 


Bo , (1) 


will exist within a superconductor when gravitational 
effects are considered, where Bọ is the external magnetic 
field and m and q are the mass and charge of a Cooper 
pair. Although gravitational effects are usually very 
weak, they may play a major role when other kinds of 
forces are nulled, as in the case of superconductivity. 
Equation (1) indicates that the residual magnetic field is 
caused by gravitational effects, since the magnetic at- 
tenuation coefficient 


2 
m'u 
B=— ~= (2) 
q'u 


is entirely determined by the ratio of the gravitomagnetic 
and magnetic permeabilities u, and u of superconduc- 
tors. Thus the essence of the problem of magnetic at- 
tenuation of superconductors reduces to the determina- 
tion of values of u and u, for superconductors. This to- 
pic forms the substance of the rest of the paper, and a 
theoretical model relating o and u is developed in order 
to predict the effects of gravitation. 


Il. CURRENT DENSITIES 


In our approach a superconductor is treated as a 
many-particle system. In the presence of applied elec- 
tromagnetic and gravitoelectromagnetic potentials A, 
®,A,,®,, the Lagrangian can be written as 


L = Sim jv7-DY9q(@-v;- A)—dm,(®, —v;- A, ) a 
j J j 


(3) 


The canonical momentum for the jth particle computed 
from the Lagrangian is related to its velocity by 


When this generalized momentum operator acts on the 
superconducting condensate wave function #7*, it yields 
the macroscopic averaged total charge- and mass-current 
densities j> and j}, respectively: 


= q jn 
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Equation (5) is the usual quantum-mechanical equation 
for a many-particle system,° where each particle has 
charge q; and mass m,; w is the superconducting order 
parameter with a phase ¢,¥=|y\e'?, which forms a su- 
perconducting condensate coherent wave such that the 
local density of superconducting electrons is given by 
n'=|wpl|*. Equation (6) is the gravitational analog of Eq. 
(5), which must apply since the charge and mass belong 
to a common carrier. For a superconductor one need 
only consider two kinds of charged particles: the Cooper 
pairs and lattice ions. 
With this information and the Maxwell equations 


5489 ©1992 The American Physical Society 


5490 


1 ; aP 3 A 
—VXVXA=j -+VXM+——e i 
Ho Jef ðt ° 3t? 

2 
ð A, 


a , (8) 
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which explicitly show all current sources, u and u, as- 
sume their free-space values. With this approach we are 
able to interpret the following two features of the macro- 
scopic equations for the vector potential A and its gravi- 
tational counterpart A, as it pertains to superconductivi- 
ty. 

(i) BCS theory shows that the superconducting elec- 
trons in the Fermi-sea ground state are expected to result 
in the formation of electron pairs with zero angular 
momentum and spin.’ Those so-called Cooper pairs 
comprise the superconducting current carriers in a super- 
conductor. We can therefore identify the macroscopic 
average of the “free” superconducting charge and mass 
currents 


iQ atta, | (9) 
and 
m m 
br OO | Acre A (10) 
Sf q q £ 
with that carried by the Cooper pairs, where the Cooper- 
pair mass m =2m,, the charge q = — 2e, and 
2 
=1 |4]? (11 
Q=-—|v ) 


is the kernal function. The Cooper pairs can move freely, 
and since their motion is driven entirely by the applied 
fields, we refer to these currents as free currents. 

(ii) The macroscopic averaged bound charge and mass 
currents, which we identify as the sources of the magneti- 
zation and gravitomagnetization, appear in Eqs. (7) and 
(8) via the quantities M and L, where M and L (Ref. 8) 
represent the macroscopic averaged magnetic dipole mo- 
ment and angular momentum densities that appear in su- 
perconductors when external fields are present. Because 
the Cooper pairs possess no angular momentum and 
hence no magnetic dipole moment because of the fact 
that the electron pairs occupy the ground state with 
equal and opposite spin, the only possible contributors to 
M and L are the lattice ions such that 


S cot OE er a . 
jem =~VXm 2im, (Y Viv WV wv Fs (12) 


Inm Y XL/2= ŽV, pyy") , (13) 


where q, and m, are the charge and mass of a lattice ion, 
and V, is only related to the coordinates of the lattice 
ions. It follows from the Le Chatelier principle? that the 
superconducting system should exert a stronger influence 
on the magnetic fields which aim to destroy it. The 
Cooper pairs reacting to the effect of a magnetic field less 


NING LI AND D. G. TORR 46 


than B, do not dissociate, but their center of gravity ac- 
quires considerable angular momentum, which gives rise 
to a superconducting condensate wave, and this induces a 
larger internal magnetic field. The interaction energy of 
this internal magnetic field with the magnetic moment of 
the lattice ions drives the lattice ions and superconduct- 
ing condensate wave function to move together vortically 
within the range of the coherent length and results in an 
induced precession of the angular momentum of the lat- 
tice ions. Consequently, a time-dependent gravitomag- 
netic fields is generated within the superconductor until 
the condition B+ (m /q )B, =0 is satisfied. 

Then the total charge and mass currents given by Eqs. 
(5) and (6) can be replaced by 


; q? m 
je =2M, |4’ V, 9—4 A ne Yl? , (14) 
M 
in =27 lV,- A+ TA, wl, (15) 
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where M,=q,2/2m, is the Bohr magneton of a lattice ion 
and y, =q,/2m, is its gyromagnetic ratio. 


IHI. MAGNETIC ATTENUATION COEFFICIENT 


Usually, the relative magnetic and gravitomagnetic 
permeabilities 


Ga 
u,=1 pee (16) 
Ger) 
and 
Tam Sin 
Hg 7l ssa (17) 
aa 


are used to represent the effect of the bound currents, 
which are dependent on the material electromagnetic 
properties as well as on the internal structure of the su- 
perconductor. Then Eqs. (7) and (8) can be simply writ- 
ten in the form 

2 


YVxVX A=pip Hote y l (18) 


2 
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a : (19) 


VxvxX A,= at en ae — Hg 0&0 


where the electric polarization P=e,y,E and the permit- 
tivity €= £€( X; +1)=e,e, are used. Since there is no such 
of the gravitational counterpart of the electric dipole mo- 
ment density P, the necessary condition required by gen- 
eral relativity, €, =€, o is satisfied, i.e., the gravitoelectric 
permittivity ¢,=1/47G is a constant, where G is 
Newton’s gravitational constant. The relationships of u 
nae Hg with their free-space values uo and Hgo, respec- 
tively, 


Jem ies 


1+ 
Giy 


H= Ho 
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(21) 


follow at once. For understanding the physical meaning 
of Eqs. (20) and (21), we may consider the following cases. 

(i) Clearly, if joy, and jm m are neglected, then Eqs. 
(20) and (21) assume their vacuum values, u=U, and 
Hg Hgo Below, we demonstrate that the macroscopic 
quantum-mechanical nature of the superconducting state 
is crucial for understanding that these bound currents 
cannot be neglected. Equation (12) indicates that each 
lattice ion will acquire momentum from the interaction 
energy of its magneton with the magnetic field of the su- 
perconducting condensate wave function in the amount 
of 


Dp, =24V,¢ , (22) 


where the density of the lattice ions is assumed to be 
equal to that of the Cooper pairs as |Ņ|?/2 and the rela- 
tion 


ia =q lY Pp, /2m, or in m= IYI P: /2 (23) 


is used. Integrating Eq. (22) around a closed path, one 
finds'° 


f p,-dl=2nh , (24) 


where the line integral of V,ġ around a closed path must 
be set to an integral multiple n of 27 to satisfy the physi- 
cal requirement that the order parameter yw be single 
valued. One finds that the vortical momentum for a sin- 
gle superconducting lattice ion must be an integral multi- 
ple of Planck’s constant h. The uncertainty principle 
provides a basis for arguing that the integral n is not 
zero. Since the product of momentum and position must 
satisfy the condition 


hi 
nies 25 
the quantum number n thus must assumes values 
n=l. (26) 


It implies that each of the lattice ions will possess uni- 
formly quantized angular momentum 


6 L,d0=2nh or L,=2nh, (27) 


where £,=rXp, and the integral element dl=r 
X(L£/L)d6. Thus the angular momentum is quantized 
in integral multiples #. We conclude that the uniform 
quantized vortical motion of the lattice ions results natu- 
rally in the fluxoid quantum of superconductivity regard- 
less of whether or not there are holes in the order param- 
eter w. 

(ii) The perfect diamagnetic property of superconduc- 
tors or the Meissner effect suggests that u,=0. It is 
equivalent to the statement that the magnetic field is zero 
inside a superconductor. However, results of experimen- 
tal measurements of the Meissner effect vary widely. 
Several theoretical models have been preoposed to ex- 
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plain the departures from the ideal Meissner effect, for 
example, the effect of trapped flux and impurities. To the 
best of our knowledge, the attenuation of magnetic fields 
by superconductors has as yet not been satisfactorily ex- 
plained. 

If the model presented in Ref. 1 is correct, then gravi- 
tational effects may be used to understand features of the 
Meissner effect, since the values of the ratio of u, and u 
can be derived theoretically by rewriting Eqs. (20) and 
(21) as 


Mgr 
H,— l 


-s 
Jm,M 


(28) 


Jm,f | | Jem 
Note that je ¢|lim,¢ and jm, mllje m- Substituting Eqs. (9), 
(10), and (11)—(13) into Eq. (28) and taking account of the 
experimental fact that u, <<1, we obtain 


Del = 8 aE, (29) 
Ho H, 


where y, =|q|/2m is the gyromagnetic ratio for a Cooper 
pair. It should be pointed out that since nothing is 
known of the phase velocity v, of a gravitational wave 
propagating within a superconductor, it is usually 
presumed to be equal to the velocity of light, c. We argue 
that the interaction of the coupled electromagnetic and 
gravitoelectromagnetic fields with the Cooper pairs in su- 
perconductors will form superconducting condensate 


waves characterized by a phase velocity v,. Since 
HoE ~ l/c? =p, o€¢.0 and pwe=k*/w’ =p, €, =U. the 
phase velocity can be predicted for the first time as 
eee eT. 9 
== 10°(m/s) , (30) 


which is 2 orders of magnitude smaller than the velocity 
of light, c. This difference arises because £, =€, o and the 
value of the gravitomagnetic permeability yw, 
=(7./Y1)Mgo given in Eq. (29) is y,/y, orders of the 
magnitude larger than its free-space value u, o~ 10° *°. 

Finally, one finds that the magnetic attenuation 
coefficient 


2 
mM Hgo Ye 1 
Bz, = | ma 
q Ho Vr Ur 


is (y./y,)(1/u,) times larger than its free-space value of 
Bo=m iby o/ q°uo~ 107} and its magnitude is deter- 
mined mainly by the magnetic permeability of the super- 
conductor. For example, if Shiffs screening factor!! 
B~10~’ and y, =3.24X10’ s~'T~! for °H, are used in 
Eq. (31), then one finds that the permeability of `H, is the 
order of u, ~ 10~ °°, which is small enough to have a very 
significant consequences for the internal magnetic field. 
Our finding is not only supported by Schiffs theoretical 
screening factor, but also experimentally supported by 
the measurement of minimum remnant magnetic fields of 
B~10~'* T,* which corresponds to attenuation of the 
Earth’s magnetic field by a factor of ~ 1077. 
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IV. o-u MODEL 


The measurement of u for superconductors is difficult, 
because its values are near zero. Thus a theoretical mod- 
el that relates u to some measurable quantity is needed. 
Below, we develop a model which functionally relates u 
with the conductivity ø. Choosing the London gauge for 
the vector potentials 


(V-A)=0, (V-A,)=0 (32) 


and introducing the Fourier components of the appropri- 
ate quantities in Eqs. (9), (18), and (19), we obtain 


je (k,@)= — Q(k) Ako) t Ako), (33) 
H! 
k? A(k,@) 
H, 
=—uQ(k) | A(k,@)+— A,(k,@)| , (34) 
k? | Zer A,(k,@) 
Leg, 
= Hg Q(T A(k,o)+— A,(k,@)| , (35) 
k: A(k,o)=0, k-A,(k,o)=0, (36) 
which yields the result 
A m Ree g 
iat 2 A (37) 
A qu H, Hgr™l 
Equating two expressions for the current, 
je(k,o)=o(k,œw)E(k, w) 
m? ‘a r 
=—Q(k) {1-8 Fert alk), 
q H H, lz l 
(38) 


we can define a complex conductivity proportional to 

Q(k) such that 

—iQ(k) ne mp, Lt, Le, 
a qu M, Hgr™l 


Before we examine further the relationship between u, 
and o, we consider the expression for the Kernel function 
given in Eq. (11). One can see that the kernal function is 
independent of k. Thus the complex conductivity [Eq. 
(39)] can be written as a function of œ only of the form 
a(w)=a,(w)—io,(@). Both the real and imaginary 
parts of o (œ) enter into the determination of the response 
of a superconductor to time-dependent electromagnetic 
and gravitoelectromagnetic fields, as shown in Fig. 1.’ 
For a given E field, the energy absorption per unit 
volume is determined by the real part, i.e., o ;(@)E”. One 
sees from Fig. 1 that at low temperature o ,(@) falls ex- 
ponentially, and for Aa/2A <1 (A=1.73%T.), o,(@) is 


o(k,o)= 1 (39) 


hw 
2A 


FIG. 1. Complex conductivity of superconductors in the ex- 
treme anomalous limit (or extreme dirty limit) at T=0. The 
rise of 0, as 1/w below the gap describes the accelerative super- 
current response (Ref. 7). 


equal to zero. Of particular importance is the propor- 
tionality of 0,(@) to 1/@, as shown in Fig. 1. One sees 
that T—0 as oa—0 and 0,(w)*1/w@—o. In fact, in the 
range of iw /2A < 1 there are no thermal excitation mech- 
anisms present, and the only process allowing absorption 
of energy is the creation of the Cooper pairs, which is 
determined by only 0,(@). This allows one to drop the 
real part and to write the conductivity as 


= Ma we gr 


q°u H, Hgr` l 


1 zt 
ee “pP 
No) ao 47 


| j (40) 


where w*,=41n‘q’/m is the plasma frequency. One may 
write Eq. (40) in the standard form for the temperature- 
dependent conductivity, !? namely, 


2 
m 
—> 7 ES eee aa 


: (41) 
m Heo Ye 


obtained by integrating o(q@) in Eq. (40) over frequency 
space, where the expression for 7 is given by! 


2 
; (42) 


| 


X 
sinhx 


=4rk,T f° <2 a%Xo)Fo) 


T 


with x =fiw/2k,T. The expansion of x /sinhx for x <1 
1S 


x l 
sinhx = 1+x?/6+x*/120+ --- 
2 4 
pg ess a init (43) 
6 120 
It yields the coupling constant 
at =2 f° £9920) F(@)[1—O(x?)] (44) 
0 O 


Substituting Eqs. (42) and (44) into Eq. (41), it follows at 
once that 
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_ p h 


_ GM Mo Yi 
4r 2mrk,TA* 


p (45) 
m Hg,0 Ye 


o(T 


If the gravitational effect is neglected, i.e., if there is com- 
plete shielding, i.e., B=O. Equation (45) then becomes 
the well-known formula for the conductivity, 


_ nq? 


m 


oo To » (46) 


where the relaxation time 


h 


Sic (47) 
27k, TA* 


To 


describes the time constant for the supercurrent to die 
away when the applied fields are moved. Its value is 
determined by the temperature T and the coupling con- 
stant A*. In copper, for example, n°=8.5X 10? m~?, at 
T =273 K, the conductivity 09>=6.4X 10’0~'!m™', and 
hence the relaxation time 7=2.68 X 10° '4s. However, at 
T =4 K, copper has a conductivity 606 times higher, so 
that the relaxation times becomes 7)=1.62X10 !! s. 
This result would appear to indicate that the copper is 
not a superconductor even at very low temperatures. 
Thus, except for the condition of low temperature, super- 
conductors must be characterized by some other unique 
properties which distinguish them from normal metals at 
low temperature which allows the near-infinite conduc- 
tivity. 

Equation (45) suggests that one of these unique proper- 
ties is the near-zero permeability. Rewriting the conduc- 
tivity approximately as 


(48) 


one sees that u,—>0, 0 < 1/u?—> œ. Estimating the order 
of magnitude of the terms in the above equation, 
mu, o/q°My~O(10-*), = -y./v, ~O(10*), and ao 
~ O(10!°), one finds that 


el 
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o~ 1072- : 
Hr 


(49) 


Substituting published measurements of the conductivi- 
ties!? (which are at least larger than 2.4X 10% Q7!m7™! 
for superconductors) into Eq. (49), one finds that the ex- 
perimental measurement of œ indicates that 
u,~O(10~7’), which is of the order of 10° larger than 
the theoretical value obtained by using Shiff screening 
factor for °H,. 

There exists, however, an inconsistency of the relaxa- 
tion time as derived from conventional superconductivity 
theory. Equation (48) also can be simply written in the 
form of the relaxation time 


m'Hgo Yi 1 


q'Ho Ye p? 


(50) 


TTo 


If t9>=%/2rk, TA* ~O(107"') and if u, ~0(1077), Eq. 
(50) gives r~O(10*) s. Although this result is of order 
10” longer than that of copper at low temperatures, it is 
still not consistent with the extremely long persistence 
time of supercurrents. Previously proposed resolutions of 
this apparent inconsistency are that at any nonzero fre- 
quency there will be an ac noise voltage reflecting the real 
part of the ac impedance of the superconductor. Then a 
superconductor is really a perfect conductor only for 
direct currents, since, as w—0, oo. We argue that 
this proposed resolution is not an experimental fact. 

We suggest that this inconsistency may be caused by 
taking integral of Eq. (40) over all frequency space. Since 
the frequency w is related to the temperature, when in- 
tegrated over the whole of frequency space, the result ac- 
tually includes the nonsuperconducting range of conduc- 
tivity as well. Therefore an additional equation is re- 
quired in order to limit the temperature effect as well as 
to isolate the magnetic permeability effect. In fact, if M 
and L are explicitly included as in Eqs. (7) and (8), Eqs. 
(33)—(35) give 


ik x —4*— (ko) — kX — 2 — ™ iko) 
m ,—1 Hg l 2q 
Alko) +" A,(k,o)= = , (51) 
m ro r 
pol a Hr Hg! 
ikXuoMlk o) — 1X go (0) FE 
ji(k,o)=—Q(k) 2: am a (52) 
r m E~ r 
ktp |) — Heo” h He, 
Hu —1 Hog H, Hyd 


The above two equations yield the needed additional 
equation 


O Mug Mp1 Hg, 


q'u Hr Hgr` 1 


Ew = 0) 


Aa (53) 


> 


Combining Eq. (53) with Eq. (39) to eliminate any ex- 
plicit dependence on œw, the major role played by the 
gravitational effect in superconductivities becomes ap- 
parent, since we find the following simple relation be- 
tween o and u, which we refer to as the o-u model: 
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Conductivity vs Permeability For Nb 


Conductivity øo (10% Siemens/m) 


Permeability u, 


FIG. 2. o-u model: conductivity as a function of permeabili- 
ty computed from Eq. (54) for several value of density of the 
Cooper pair of the order from 10” m~? (bottom line) to 10 
m`? (top line) for niobium. 


Fa M, €00 
g= ———— 
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mMM, Hg o 1 


i 
qho My 


+1 z 


(54) 


One sees that the conductivity is only determined by the 
magnetic permeability or vice versa. The conductivity 
calculated from Eq. (54) as a function of permeability and 
the density of the Cooper pairs for Nb is illustrated by 
the numerically computed logarithmic curves in Fig. 2. 
It can be seen from Fig. 2 that the conductivity increases 
exponentially with continuously decreasing permeability. 
This characteristic dependence can be used to obtain the 
values of u experimentally. For example, Fig. 2 shows 
that the experimentally measured values of conductivities 
larger than the order of 10%% are consistent with per- 
meabilities smaller than 107% over densities ranging 


Relaxation Time vs Permeability For Nb 


Relaxation Time T S 
Oo 
on 


1073" 10730 


10778 10728 10724 
Permeability u} 
FIG. 3. Relaxation time dependence of permeability for vari- 


ous densities of the order from 1072 m~? (bottom line) to 10? 
m`? (top line) for niobium calculated based on Eq. (55). 
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from 10” to 10? m~? for a Nb superconductor. 

The implications of the o-u model are made clear by 
returning to consideration of the relaxation time intro- 
duced in Eq. (46). Equation (54) can be written in the al- 
ternative form 


1/2 
mM, Eo 


mQ 


_ il 
ere 


l, 


MM Ugo 1 


q'Ho H; 
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By varying u, from 107°? to 107%, the numerically com- 
puted curves of relaxation time as a function of permea- 
bility based on Eq. (55) are those illustrated in Fig. 3. It 
is shown that the relaxation time increases exponentially, 
and for the typical case of u,<10~7°, the supercurrent 
can persist for at least years. This extremely long lifetime 
for the supercurrent is indeed consistent with the experi- 
mental observations of the persistence of superconduc- 
tivity. 


V. CONCLUSIONS 


Expressing the conductivity in Eq. (54) as a function of 
the magnetic attenuation coefficient as shown in Fig. 4, 
one sees that for typical observed conductivities larger 
than the order of 10%, the corresponding magnetic at- 
tenuation coefficient is larger than the order of 107%. 
Comparing this experimental limit with the extreme 
theoretical approximation of Schiffs result, a possible 
range for the magnetic attenuation coefficient can be pre- 
dicted as 


Conductivity vs Attenuation For Nb 


Conductivity ø (10% Siemens/m) 


Attenuation £ 


FIG. 4. Conductivity as a function of the magnetic attenua- 
tion coefficient calculated from Eq. (55) for different densities of 
the order from 10? m`? (bottom line) to 10? m~? (top line) for 
niobium. 
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10°7>B>10~ 2°, (56) 


The result suggests that the experimentally observed very 
poor Meissner effect* may not be simply caused by im- 
purities alone, but may be the result of a more fundamen- 
tal physical mechanism, namely, gravitational effects. 
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